Entropy at the Outskirts of Galaxy Clusters as Implications for 
Cosmological Cosmic-Ray Acceleration 

Yutaka Fujita 

Department of Earth and Space Science, Graduate School of Science, Osaka University, 

Toyonaka, Osaka 560-0043, Japan 

f uj itaOvega. ess . sci . osaka-u. ac . jp 

and 

Yutaka Ohira and Ryo Yamazaki 

Department of Physics and Mathematics, Aoyama Gakuin University, Fuchinobe, 
Chuou-ku, Sagamihara 252-5258, Japan 

ABSTRACT 

Recently, gas entropy at the outskirts of galaxy clusters has attracted much 
attention. We propose that the entropy profiles could be used to study cosmic- 
ray (CR) acceleration around the clusters. If the CRs are effectively accelerated 
at the formation of clusters, the kinetic energy of infalling gas is consumed by 
the acceleration and the gas entropy should decrease. As a result, the entropy 
profiles become fiat at the outskirts. If the acceleration is not efficient, the 
entropy should continue to increase outwards. By comparing model predictions 
with X-ray observations with Suzaku, which show fiat entropy profiles, we find 
that the CRs have carried < 7% of the kinetic energy of the gas away from the 
clusters. Moreover, the CR pressure at the outskirts can be < 40% of the total 
pressure. On the other hand, if the entropy profiles are not fiat at the outskirts 
as indicated by combined Plank and i? 05*71 T observations, the carried energy and 
the CR pressure should be much smaller than the above estimations. 

Subject headings: cosmic rays — galaxies: clusters: general — galaxies: clusters: 
intracluster medium — large-scale structure of universe 



1. Introduction 



The existence of cosmic-rays (CRs) in clusters of galaxies is clearly indicated by radio 

synchrotron emission (iFeretti et al.l 120121). The CRs may be accelerated at shocks or turbu- 

lence in the X-ray hot gas (e.g. Fujita fc Sarazin 2001; Brunetti et al.|[2001 : Fujita. Takizawa. &: Sarazin 
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20031 : iKang fc Joneal2005l : iPfrommer et al.ll2007l ). and they may play i mportant roles in the 



evolution o f clusters. For example, they may heat the cores of clusters (iFujita fc Ohirall2011 



2OI2I . l2013l ). Since the background gas is hot, their acceleration efficiency in clusters could 
be very different from that of objects in the Galaxy (e.g. supernova remnants; SNRs). How- 
ever, the lack of firm observational evidence of CRs with the exception of radio observations 
prevent us from understanding their acceleration. Thus, other observational approaches that 
can clarify CR acceleration in clusters are strongly desired. 

X-ray observations with Suzaku have revealed the entropy of the hot gas of clusters at 
their outskirts. Here, the entropy is defined as i^' = fc^T/n'^s"^, where ks is the Boltzmann 
constant, 7^ (= 5/3) is the adiabatic index for the gas, and T and n are the gas temperature 
and the number density, respectively. Numerical simulations, in which non-gravitational 
heating (e.g. energy ejection from galaxies) and radiative cooling are not included, predicted 
that th e entropy should increase with r as oc r^'^, where r is the distance from the cluster 
center ( IVoit. Kay, fc Bryan! l2005l : iBurns. Skillman. fc O'Sheall2010l ). On the oth er hand, the 



Suzaku observations have shown that the entropy p r ofiles fiatten at r > 0.5 r9nn (IBautz et al. 



Kawaharada et al. 


2010; 


Hos 


uno et al. 


2010; 


Akamatsu et al. 


2011; 


Simionescu et al. 


Sato et al. 


2012; 


Walker et al. 


2012a: 


Ichikawa et al. 


2013 


) (see Figure 13 of 


Sato et al. 



)12l ). where r2oo (~ 2 Mpc) is the radius within which the mean mass density is 200 times 
the critical density for a flat universe and is often considered to be the typical radius of 
clusters. 

Several ideas have been proposed to solve this discrepancy. Gas clumping could lead 
to the overe stimation of the ga s density, which in turn could lead to the underestimation 



of entropy (INagai fc Laul I2OIII ). However, gas clumping is effective mainly at r 



> 



^200, 



because the clumps tend to be destroyed within the cluster (r < r2oo)- Non-equilibrium 
between ions and electrons may also decrease entropy. This is because heavier ions have 
most of the kinetic energy of gas and they are first thermally heated at shocks. If it takes 
a long time for the kinetic energy of the ions to be transferred to lighter electrons, the 
observe d temperature of the electrons could be low. This may decrease the observed entropy 
of gas (IHoshino et al.l I2OIOI ; lAkamatsu et al.l 1201 ll). However, nurnerical simulations have 



shown that the degree of non-equilibrium is small (IWong fc Sarazin l2009h. although it ma y 



depend on the dynamical state of clusters (lAkahori fc Yoshikawa 



2008 



Rudd fc Nagaill2009h . 



Moreover, from observations of the Sunyaev-Zel'dovich effect, Planckhas shown that the gas 
pressure at the outskirts of clusters does not significantly drop. Th is suggests that the 
ions and electrons are in equilibrium (jPlanck Collaboration et al.ll2013l ). For some clusters, 
accretion of matter toward them may decrease at low redshifts. Although this makes the 
entropy profiles fiatter at their outskirts, it does not seem to be effe ctive enough to be 
consistent with the observed systematic fiattering (jCavaliere et al.l 1201 ll ). 
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In this Letter, we propose that the entropy profiles may refiect CR acceleration in clus- 
ters and could be used as a tool to study it. We focus on the acceleration at 'accretion shocks' 
around the clusters. These shocks are generated through gas inf all toward the clusters as they 
grow by accreting dark matter and gas fr om the outside (e.g. lEvrard. Metzler. &: Navarro 



19961 : iRyu et al.ll2003l : ISkillman et al.ll2008l ). If the CR acceleration at the shocks is effective, 



the CRs sh ould consuni e much of the kinetic energy of the gas and should affect the entropy 
of the gas (lBykovll2005l ). Thus, the flattened entropy proflles may be the result of the past 
CR acceleration around the clusters. 



2. Models 



It is widely beheved that if C Rs are effectively accelerated at a shock, their pressure 
chang es the structure of the shock (jPrury fc Voelklll98ll : lDrurylll983l : iMalkov fc O'C Drury 
200ll ). Before we investigate the entropy proflles of clusters, we consider the jump conditions 
at the shock. We treat gas and CRs as two fluids hereafter. We deflne three distinct regions 
with a shock: (0) the far upstream region where the influence of accelerated CRs can be 
ignored, (1) a region just upstream of the shock where the presence of the CRs affects the 
shock structure, and (2) a region downstream from the shock. We indicate each region by 
lower indices. For example, the gas density and velocity in the far upstream region are po and 
f 0, respectively. The spatial scale of the region (1) is much smaller than the size of a cluster, 
and it cannot be resolved by current X-ray telescopes. The gas density and temperature 
could respectively change from p2 and T2 far downstream of the shock as the cluster grows. 

The CR acceleration consum es gas energy. We discuss the accompanied entropy change 
at a shock based on a model of IVink et al.l ( I2OIOI ). We represent the energy flux that is 
carried away by CRs diffusing away far upstream by F^r- The flux is normalized to the 
kinetic energy flux of the shock: 

cr / ^ N 

The fraction of the downstream CR pressure is 



w 



2,cr 



2,cr 



P2.th + F. 



(2) 



2,cr 



where 'cr' and 'th' represent the CR and the thermal components, respectively. The Mach 
number of the shock deflned at the region far upstream is 



Mo 



7; ^0 



(3) 
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The compression ratios across the different regions are 



Pi P2 
Xi = —, X2 = — , Xl2 
Po Pi 



P2 

X1X2 = — ■ 
Po 



(4) 



Assuming that the thermal component evolves adiabatically in the region (1), the pressure 
is given by Pi^th = -PoX?- Thus, the total compression ratio is given by 



X12 



ha + ^)Mlxi 



79 



(7, - l)Mlxi 



(79+1) 



+ 2 



(5) 



The CR pressure fraction is given by 

[l-X?)+laMl{l-l/xi) 



w 



1 + 73^42(1 -I/X12 
and the escaping energy flux is given by 

2Gq 2G2 2(j2 



(6) 



where 



1 + 



I9 



lgMlXl2 X12 



1 - 



X12 



2 ' 
XI2 



Go = w- 



7c 



+ {l-w)- 



19 



(7) 



(8) 



75-1 7cr - 1 79-1 

and 7cr is the adiabatic index for the CR component. Equations. ©-([H]) show that eesc and 
w are represented by xi for a given AAq. The entropy difference across the shock is written 
as 

'K2 



In this equation, 



K2 



:i-w] 



79 
A12 



which means that AS* is also the function of Mq and xi ( jVink et al.ll2010l ). 



X12). 


5 


fVink et al. 


2010) 



(9) 
(10) 



I t would be natural to assume that CRs are a ccelerated at the accretion shocks (IFujita fc Sarazin 



2OOII : iKang fc Jones! I2OO5I : iPfrommer et al.ll2007l ). although clear evidence has not been dis- 



covered. The CRs would carry the gas energy away and lower the entropy. We compare the 
entropy profiles predicted by numerical simulations that do not include CR acceleration with 
those obtained through X-ray observations. Contrary to density and temperature, entropy 
is conserved if there is no heating and cooling, which is an advantage to discuss entropy. The 
entropy profile predicted by the numerical simulations is 



1.32/^200(^^^200) 



1.1 



(11) 
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where ^^^200 ~ 362 keV cm^(Tx/keV) and Tx is the gas temperature (IVoit et al.ll2005l ). If we 
adopt the entropy profiles obtained with Suzaku, they can be fitted as 



(12) 



(jWalker et aLll2012bl ). At r ~ 0.3 r2oo, the entropy profiles do not become fiat and they 
are not much afi^ected by the activities of central active galactic nuclei (AGNs). Thus, we 
assume that 

-f^obs(0.3r2oo) = ^sim(0.3r2oo) , (13) 
which gives the normalization of Equation (fT2!) . Al though there is a d ebate about an ap- 
propriate position of the normalization (r ~ 0.3 r2oo; lEckert et al.ll2013l ). it would not affect 
the following results when the observed entropy profiles significantly deviate from those pre- 
dicted by the theoretical simulations. Figure [T] shows -ft"sim('") and Kohs{r) for Tx = 8 and 
4 keV. Although Kohs slightly decreases at r > 0.7r2oo, the curve should be regarded as 'fiat' 
for r > 0.5 r2oo considering the errors of the observations. We not e that while non-thermal 
press ure owing to turbulence in gas may affect the entropy profile (ILau. Kravtsov. &: Nagai 
20091 ). it should be included in -R'sim obtained by high- resolution numerical simulations that 
can follow turbulence. Thus, the difference between Ksim and i^Tobs is made by something 
other than the non-thermal pressure. 

We assume that a cluster gradually grows from the inside to the outside as it accretes 
matter from the outside. The entropy of the gas that falls into the cluster may not be 
zero, because it could have been heated through activities of galaxies and AGNs in the field 
(Kaiser 1991 : Cavaliere. Menci. fc Tozzi 1997 : Ponman. Cannon, fc Navarrol 1999[). In fact , 
widely distributed metals outside clusters may indicate such activities (IFujita et al.ll2008l ). 
The entropy of the gas befo re passing accretion shocks is estimated to be Kout ~ 100- 
400 keV cm2 Jvoit et al.ll2003l ). and we adopt K^ut = Kq = 200 keV cm^. 



We ignore heating and cooling of the gas except for the heating at accretion shocks 
around the cluster (r > r2oo)- The motion of the gas is fundamentally controlled by the 
gravity from dark matter, and thus the infall velocity of the gas toward the cluster and the 
Mach number of the accretion shocks (Aio) do not significantly depend on the presence or 
absence of CR acceleration at the shocks, while the postshock gas (region 2) is significantly 
affected by it. Therefore, the Mach number of the shock Aio = Ai{r), through which the 
gas located at the radius r at present had passed, can be estimated using Ksim{T) and Kout- 
Since i^'sim('") is the entropy when the CR acceleration is ignored, the relation among Ai, 
i^'sim and Kout should be described by the Rankine-Hugoniot relation that does not include 
the effects of CRs (xi = 1): 



i^sim(r) _ 2-fgM{ry - (7, - 1) 



K 



out 



7, + ! 



(7g-l)A^(r)^ + 2 
(7, + l)M(r)2 



79 



(14) 
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Figure [2] shows the profiles of the Mach number A4{r) derived from Equation f|T^ . The 
Mach number is the increasing function of r and reaches ~ 10 at r ~ r2oo- We emphasize 
that entropy is conserved for a given gas element, and that the gas observed at r at present 
was not necessarily there when it passed a shock. Therefore, the radius of the shock could be 
larger than r2oo «^ the time when the gas passed the shock. In other words, the position of 
the shock is not specified in our model. 

Assuming that CR acceleration makes K^hs smaller than -ft'sim at the outskirts of clus- 
ters, the entropy difference at the shocks for the observed clusters is obtained by using 
Equations @j and ([10]) with Mo = M{r): 

Kout 

Since we already know _K'obs('") (Equations [I2] and [I3]) and Ai{r) (Figure [21), we can derive 
Xi(r) from Equation f[T5]) . 



AS[Mir),Xiir)] = hB\n 



(15) 



Results and Discussion 



FromA^o = A^(r) and Xi(r), we obtain eesc('") = ^esc[-M{r),Xi{r)] andw(r) = w[Ai{r),Xi{r)] 
using Equations ([S])-®, which are shown in Figure [HI In the calculations, we assumed that 
7cr = 4/3. Figure [3] shows that egsc and w increase outwards and reach e p^^ ^ 0.07 and w y 0.4 
at r ~ r2oo- These values are often estimated for the Galactic SNRs (IVink et al.ll2010l ). al- 



though their Mach numbers are typically an order of magnitude larger than those of the 
cluster shocks considered in this study. Nevertheless it is interesting t hat clusters have si mi- 
lar values of eesc and w. Since other factors such as cl umping of the gas jNaeai fc Laulboilh or 
a decrease of matter accretion toward the clusters f Cavaliere. Lapi. fc Fusco-FemianJ2011 ) 
may partially contribute to the decrease of the entropy, the actual values of eesc and w could 
be somewhat smaller than the above values. 

We implicitly assumed that the accretion shock is a single shock. However, the accreted 
gas may have passed multiple shocks. In that case, the Mach number of the inner shocks 
may be smaller than in Figure [2l because of the heating at the outer shocks. Thus, 
our model could not be applied to dynamically active clusters that have complicated shock 
structures. However, the CRs may be reaccelerated at mult iple shocks and the acceleration 
efficiency may be large in spite of the low Mach numbers (iKang fc Ryull201ll ). Therefore, 
CR acceleration at multiple shocks with small Mach numbers may be qualitatively similar to 
that at a shock with a large Mach number. At least we can say that the difference between 
i^obs and Ksim indicate that a significant fraction of the gas energy must have been converted 
to CRs even in the case of multiple shocks. 
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Some cosmological numerical simulations do not show the drastic temp erature decrease 
or density increase at cluster outskirts in spite of CR acceleration (e.g. iPfrommer et al. 



20071 : IVazza et al.ll2012l ). Although we could not identify the cause, it may be because their 
assumed acceleration efficiency is smaller than those we obtained. In fact, i f a constant high 



efficie ncy of 50% is adopted, the temperature significantly drops (Figure 3 of iPfrommer et al. 



20071 ). Moreover, the efficiency in the central region of a cluster may need to be much lower 
than that at the outskirt in order to reproduce the flat entropy proflle. In our model, 
the CR acceleration is prohibited in the central region because of the low Mach number 
attributed to the preheated gas (-K'out > 0). This leads to a low density of CR protons and 
less interaction of CR protons with gas protons in the central region. This may be favorable 



for th e non-detection of gamma-rays from the central regions of clusters ( lAckermann et al. 



20101), because gamma-rays are produced through proton-proton interaction. 



Our model predicts that entropy decreases in the region where CR acceleration is ef- 
fective. The CR pressure is w > 0.1, which is the typical value for the Galactic SNRs 
( IVink et al.ll2010l ). at 0.5 r 200 ^ ^' ^200 (Figure [3]d). I n this region, synch rotron emission 
from electrons ('radio relics') has often been discovered ( iFeretti et al.ll2012l ). Although the 



electrons seem to be accelerated at shocks created during cluster mergers, the Mach num bers 
are relatively small {M. ~ 2-4) (Ivan Weeren et al.ll2011t lAkamatsu fc Kawaharall2013l ). At 
these small Mach numbers, CR acceleration may be difficult. However, if there is a pre- 
existing CR population, CR reacceleration at the shocks could increase their energy high 
enough to emit the synchrotron radiation ( iKang fc Ryul I2OIII ). The entropy decrease at 
cluster outskirts may indicate the existence of a pre-existing CR population. Moreover, the 
CRs r nay amplify magii etic fields ( iLucek fc Belli I2OOOI : iBelll 120041 : iFujita. Ohira. fc Takahara 
2OIOI : iFujita et al.ll201ll ). Thus, magnetic fields amplified by the CRs may be found at the 
cluster outskirts in the future. The CRs accelerated at the accretion shocks may contain 
ultra-high-energy ones ( llnoue. Aharonian. fc Sugiyamal 120051 ). 



Note that the flat entropy proflles have not been conflrmed in the combined analysis 
of the pressure pro flles obtained with Plank and the density proflles obtained with ROSAT 
( lEckert et al.ll2013[ ). In that study, the entropy continues to increase outward. If this is 
the case, our model indicates that the CR acceleration efficiency is much smaller than those 
estimated above. This could constrain the aforementioned reacceleration model for the radio 
relics. This also shows that the future confirmation of the entropy profiles is important in 
terms of CR acceleration in clusters. At present, the model and observational uncertainties 
suggest that the estimated values of eesc ~ 0.07 and w ~ 0.4 at r ~ r2oo at the beginning of 
this section should be regarded as upper limits. 
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Fig. 1. — Profiles of -ft'sim (thin lines) and K^hs (thick lines) for Tx = 8 keV (solid lines) and 
Tjsf = 4 keV (dashed lines) 
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Fig. 3. — Profiles of (a) eesc and (b) w for Tx — S keV (sofid lines) and Tx — 4: keV (dashed 
fines) . 



